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ABSTRACT 

Grid stars and reference stars provide the fundamental global and local 
astrometric reference frames for observations by the Space Interferometry 
Mission. They must therefore be astrometrically stable at the ~ 1 /zas 
level. I present simple formulae in closed form to estimate contamination 
of these frames by motions due to stellar companions that go undetected in 
a radial-velocity (RV) survey of specified precision. The contamination rate 
depends almost entirely on the binary-period distribution function and not on 
the details of the mass or eccentricity distributions. Screening by a modest RV 
survey (a rv = 60 ms -1 ) can reduce the fraction of grid stars having detectable 
unmodeled accelerations to <C 1%. Reference- star selection promises to be much 
more challenging, partly because the requirements of astrometric stability are 
more severe and partly because the required density of sources is ~ 100 times 
higher, which latter implies that less satisfactory candidates will have to be 
accepted. The tools presented here can help design a reference-star selection 
strategy, but a full treatment of the problem will require better knowledge of 
their companions in the planetary mass range. 



Subject headings: astrometry - methods: statistical - planetary systems 
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1. Introduction 

By their nature, astrometric observations measure the angular separation of one star 
relative to another, either absolutely or projected on some axis. Hence, if attention is 
restricted only to these two stars, one cannot determine which, if either, is stationary, and 
which is moving. 

There are basically three approaches to resolving this ambiguity. The traditional 
method is to find reference stars that, based on their photometric properties, are believed to 
be so far away that their astrometric motions are small and can be estimated analytically, 
at least statistically. The HIPPARCOS mission pioneered a second, radically different 
approach. HIPPARCOS measured the angles between of all pairs of stars separated by 
about a radian. Because each star's position was measured against many others (which 
lay in sections of the sky where the parallactic motion is substantially different), it was 
then possible to determine each star's motion individually (up to a global rotation) from a 
global solution of these measurements. The global rotation was then measured relative to 
the inertial quasar frame by making use of radio stars whose positions were known in both 
systems. The FAME and GAIA missions also plan to employ this method except that the 
quasar tie-in will be done directly in the optical. The Space Intereferometry Mission (SIM) 
will use a third approach, which is a hybrid of the previous two. All SIM astrometry will 
be done relative to two special classes of stars, "grid stars" and "reference stars" . 

The grid will be composed of 1000-3000 astrometrically stable stars spread over 
the whole sky. By repeatedly measuring the relative separations of these stars over the 
T m = 5yr span of the mission to ~ 10 //as precision, their absolute parallaxes and relative 
positions and proper motions can be determined to ^ 3 //as and ^ 2 /ias, and ^ 2 /tasyr -1 , 
respectively. The resulting SIM frame can be pinned to the extragalactic inertial frame by 
observing a few quasars (NASA 1999). Absolute parallaxes, proper motions, and positions 
of all other objects can then be determined by measuring their positions relative to the 
grid. For this to work, grid stars must be astrometrically stable on the scale of a few /ias: if 
too many of them prove unstable and "drop out" of the grid, the global grid solution could 
be undermined. D. Fischer (1998 private communication) originally argued that metal-poor 
K giants at D ^ 1 kpc would make the only truly suitable grid stars because their great 
distance would reduce astrometric perturbations due to unseen companions. This view has 
now come to prevail. S. Majewski and collaborators are conducting a photometric search for 
grid candidates (Patterson et al. 1999). These will ultimately be screened for companions 
using a radial velocity (RV) survey. 

The key problem is to determine in advance what intensity of screening is required to 
produce a grid that will either have very few drop-outs, or that is structured to be able to 
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"paper over" whatever dropouts it does suffer. This problem has to date been addressed 
primarily by Monte Carlo simulations (Frink et al. 200 f; Jacobs 2000; Peterson 2001). 

Planet searches with SIM will be carried out by measuring target-star positions relative 
to nearby 1°) astrometrically stable (< 1 /xas) reference stars. The requirements for 
reference stars are qualitatively similar to those of grid stars, and consequently I will treat 
the two in parallel. Quantitatively, however, the reference stars are much more demanding. 
First, the individual measurement error relative to reference stars is 1 /ias, an order of 
magnitude more precise than the individual grid measurements. Hence the stars must be 
correspondingly more astrometrically stable. Second, reference stars must be brighter to 
allow for higher precision measurements with a comparable exposure time: Grid stars can 
be V ~ 12, but reference stars should be 1 or preferably 2 mag brighter. In principle, 
closer (and hence brighter) stars could be chosen, but, as mentioned above, closer stars are 
astrometrically less stable. Third, reference stars are required at higher density. I will argue 
in § 5.2 that of order 8 reference star candidates must be found within ~ 1 deg 2 of the target 
star, whereas the density of grid stars is <^ 10~ 1 deg -2 . Thus, the selection of reference stars 
will likely include objects that are neither metal-poor nor exceptionally luminous, since the 
density of such rare stars is simply too low. Finally, as I show in § 5.1, grid stars that 
are accelerating approximately uniformly can perform their grid functions perfectly well. 
This is not true of reference stars: if their accelerations are at all measurable then they 
degrade the sensitivity to outlying planets that are themselves detectable only due to the 
uniform acceleration that they induce on the targets. Thus, reference-star requirements for 
astrometric stability are even more demanding compared to grid stars than the relative 
precisions of measurement would seem to imply. Although the discussion here will be 
focused on SIM, I note that similar considerations apply to the choice of reference stars for 
planet searches using the Keck Interferometer. 

Systematic study of the reference- star problem is less advanced than that of the grid. 
Marcy (2001) and Shao (2001) discuss methods of selecting reference stars and give initial 
estimates of contamination levels. 

Here I present an analytic method for determining contamination of grid stars and 
reference stars by undetected stellar companions. This approach is complementary to 
that of the Monte Carlo studies mentioned above: Monte Carlo methods can model 
arbitrarily complex distribution functions in arbitrary detail, and so can potentially capture 
everything that is known about a problem. In the present case, this intrinsic complexity 
seems formidable because binary orbits are quite varied and are described by seven 
parameters. However, the analytic approach has its own advantages. Primarily it allows 
one to see explicitly how the results depend on the assumptions and on the precision of 
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the measurements. Thus, analytic modeling can be particularly valuable as a guide while 
selection criteria and methodolgy are evolving. 

Although the underlying parameter space is quite complicated, I show that the 
contamination rates can be quickly evaluated using a few simple formulae. In their simplest 
version, these formulae depend only on the period distribution of the companions and on 
their mean mass. This version assumes circular orbits, but I also evaluate explicitly the 
correction due to eccentricities, which is small. 

For the most part, only stellar companions are considered in this paper. Brown dwarfs 
are not a practical concern because of the "brown dwarf desert", the observed lack of 
brown-dwarf companions to G dwarfs (Halbwachs et al. 2000; Marcy & Butler 2000), the 
progenitors of K giants. Planets are not a major concern for the grid because metal-poor 
giants are not expected to have many, and because the demand for astrometric stability is 
not so severe. As I briefly discuss in § 5.2, planets are a significant concern for the reference 
stars. An analytic treatment of planet contamination is feasible. It would be broadly 
analogous to the one given here for stellar contamination. However, the data on planetary 
companions do not at present suffice to use these methods to make reliable predictions of 
contamination rates. I therefore make reference to such contamination only for purposes of 
illustration and defer presentation of a detailed treatment. 

In § 2, I recapitulate some well-known statistical results that are required for the 
analysis. In § 3, I analyze screening of a sample using an RV survey. This analysis 
applies equally to grid stars and reference stars. In the next two sections on astrometric 
contamination (§ 4) and implications (§ 5), I treat these two classes in parallel. In § 6, I test 
the analytic formulae derived here by "predicting" the results of a Monte Carlo simulation 
by Frink et al. (2001). I find excellent agreement. Finally, in § 7, I summarize the main 
formulae derived in the paper. 



Let F(t) = J27=i a ifi(t) be a linear combination of n trial functions fi(t), each with 
coefficient a^. If the parameters aj are fit to measurements at times tk, (k = 1, . . . , N), with 
errors o^, then the covariance matrix = cov(cij, aj) is given by (Press et al. 1992) 



2. 



Statistics Prelude 



N 



C = b 



-1 




k=i 



- 5 - 



For the special case of a two-parameter fit F(t) = a± + a 2 t, where the errors are all equal, 
cTfc = a, equation ([I]) implies that the variance of the slope is given by, 



a 2 



W(a2) " 022 = iVwW ' (2) 

where var (t) = (tl) - (t k ) 2 . Finally, if there are a large number of measurements uniformly 
spaced over an interval [— T/2, T/2], also with equal errors but now allowing for an arbitrary 
number of parameters, then equation (|I]) becomes 



N r T / 2 



For a polynomial, F(t) = J27=i V(* ~ these components are, 

N (T/2) i+j ~ 2 

a 2 (z + 3 - l)(«-l)!tf - I)' 

and 6y = otherwise. Hence, even for a cubic (n = 4), 6 decomposes into two (2 x 2) 
matrices that are easily inverted by hand. 



3. Radial Velocity Screening 

I first ask the question: what fraction of an initial sample of candidate grid stars or 
reference stars will be rejected by a RV survey of iV measurements equally spaced over a 
time T rv , each with a precision <j Tv . I assume that the candidates have mass M = 1 M , and 
focus initially on stellar companions of mass m, i.e., 0.1 M & < m < 1M . Note that for 
parameters of interest, T rv ~ 5yr, a rv ~ 60 ms -1 , essentially all companions with periods 
P <^ 2T rv will easily be detected. For example, the velocity semi-amplitude induced by an 
m = 0.1 M & companion in a P = 2T rv = 10 yr orbit is 1.5 km s -1 , so that only extremely 
face-on orbits, sini < 0.08 would escape detection. Since these constitute about 0.3% of all 
orientations (and so ~ 0.02% of all candidates), they will be of no interest here. 

Instead, I will work in the limit of uniform accelerations, the magnitude of whose radial 
component is given by, 

Gin 63m/s / m \( r \~ 2 

a r = — — cost/ = : — cost/ (5) 

r 2 1 1 yr V0.3 M Q J V30 AU J 11 K ' 

where r is the instantaneous separation of the companion and 9 is the angle between the 
separation vector and the line of sight. 
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On the other hand, from equation (Q), [and using J2k=i^ = N(N + l)/2, 
Ek=i k2 = N ( N + + l)/6 var (*0 = - 1)/12], the RV survey can measure such 

accelerations with precision 



12(iV-l)(7 rv 16m/s/ cr 



\ iV(iV + 1) T rv ~~ yr V60m/s/ V5yr 



T 



-1 r 



iV(JV + l)/(iV-l) 



1-1/2 



6.67 



(6) 



Here <r rv includes contributions from both measurement error and intrinsic velocity 
variability of the stellar atmosphere, and the expression is normalized to N = 4 
observations. I will assume that candidates will be rejected if the observed accelerations 
exceed a threshold a trir = 2a a : if the threshold were set much lower, too many candidates 
would be eliminated by statistical noise. For an ensemble of systems, | cos#| is uniformly 
distributed over [0,1]. Hence, the fraction f p that will pass the RV cut is 



f P (r, m) 
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An ensemble of companions with semi-major axis A and eccentricity e, will be found equally 
distributed in all phases of their orbits. Hence the fraction passing the RV screening is 



f p (e, A, m) 



Gm 



< Othr 



f p (e,A,m) 



l + ^e' 



[A(l-eW 
A 2 a thr ( Gm 



> Othr 



(9) 



2"/ Gm \[A{l + e)] 2 

In the intervening range, f p (e,A,m) can be written in closed form, but the expression is 
not illuminating. For the moment, I use the approximation of circular orbits. I will address 
the question of eccentric orbits in § 3.2. 



3.1. Circular Orbits 

The main progenitors of K giants are G dwarfs, whose companion distribution has 
been studied by Duquennoy & Mayor (1991, DM91). They find that about 60% of G dwarfs 
have companions, and that the log periods of these are roughly Gaussian distributed with 
(log(P/days)) = 4.85 and {var[log(P/days)]} 1 / 2 = 2.3. The companion mass distribution 
is roughly flat over the interval 0.1 M Q < m < 1M & . By comparison with this broad 
distribution over roughly 7 decades in semi-major axis, equations fl8|) and @ show that 
f p rises from nearly to 1 in less than a decade. Therefore, it is appropriate to evaluate 
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f p (m) over all A as being equal to the fraction of candidates with A > A(m) where 
f p [A(m),m] = 




30 AU — -7—^ — I . (10) 

VO.3M y V32m/s/yr; v ; 

For masses of m = 0.1 M & and 1 M Q , these values correspond to periods of P = 10 4 4 days 
and P = 10 5 ' days. That is, the entire mass range corresponds to only about half a decade 
of the DM91 binary distribution function. Hence, to evaluate the rejected fraction f re j,rv 
over the entire mass range, it is appropriate simply to adopt the period corresponding to the 
average mass, (m) ~ 0.3 M , which is equivalent to assuming that the DM91 cumulative 
distribution function is a straight line within a bin. The error so induced is less than 1%, 
much less than the Poisson noise in the measurement of the distribution function itself. 
This yields, 

' P * ,df b D2 _ 4vr 2 [A((m))] 3 



/rej,rv — / dP , P^ _ , (11) 

Jo dP G(M + (m)J 

where dfb/dP is the binary-period distribution function. For the fiducial parameters I have 
been considering, P* = 10 4 ' 7 days. Thus, substituting the period distribution of DM91 into 
equation QTTj) , and using these fiducial parameters, / re j, r v = 29% of the candidates would be 
rejected. In the scenario I have laid out here, another 5% of the remaining 71% would be 
rejected because statistical fluctuations would cause the stars to appear to accelerate at the 
2 a level even when there was no real acceleration. 



3.2. Eccentric Orbits 

To include eccentric orbits exactly would be complicated because the transition from 
equation (|]) to (|) is complicated. However, the order of the effect can be assessed by noting 
that in both these limiting regimes, the change is accounted for by a thr — ► [1 + (3/2)e 2 ]a t h r , 
and therefore making the following substitution in equation (|1C)|): 

3 

athr -> geflthr, g e = l + -(e 2 ). (12) 
Then, following through the remaining logic in § 3.1, this leads to a change in the estimate 

of /rej,rv 

y logffe (df b /dlogP)\ P . , . 

dhgP Pt ~ log(7/4) 0.11 ' { } 

where the evaluation has again been made using /& from DM91. Because this change is 
so small, I will generally assume circular orbits, but will include a 2% adjustment for 
eccentricity. 



a p 3 df b 

A/ r ej,rv = — rlOgflfe 
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4. Fraction of Astrometric Accelerators 
4.1. Grid Stars 

Of the candidates that pass the RV screening and go on to become grid stars, what 
fraction will have accelerations that are detectable astrometrically? I choose a 3a criterion 
to avoid excessive rejection due to noise. The astrometric acceleration a is given by 

Gm . . 13uas/ m \ ( r \~ 2 / D \ -1 . , 

a= — -sm6> = — 4— smO (14) 

Dr 2 yr 2 V0.3 M G J V30 AU J Ukpc/ V ; 

where D is the distance to the grid star and all the remaining quantities are the same as 
in equation Suppose that the astrometric data in one direction are fit to the form 
ip(t) = tpo + fi t + (l/2)at 2 + Kilcos[27r(t — t )/yr], where an d Hq are the position and 
proper motion at mid-mission, II is the parallax, and k and t give the phase and projection 
factor of the parallax ellipse. Since the grid will be surveyed many (e.g., iV ~ 23) times, the 
covariances of these four parameters can be evaluated using equation @. First I note since 
the parallax term is cyclic, whereas all the others are secular, 6 i4 / y/bub^ 1. Thus, the 
parallax terms decouple and so can be ignored for present purposes. The remaining terms 
are then given by equation (H), and the resulting covariance matrix is therefore, 
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where T m is the mission duration and cr — o/yiV is the position error for the case where 
the data are fit to uniform motion. Hence, the acceleration error is 



Or 



2 /xas / cr \ / T, 



V2 uas/ Voyr/ 



yr 2 V2 fi&s/ V5 yr, 

As stated above, the threshold of detection is a t hr = 3a a . We can now define a figure of 
merit K for the relative sensitivities of the astrometric and RV surveys 

K = ^ = L13 ( ^ ) ( ath : yi^-Y 1 . (17) 
Da thr V32 m/s/yr/ V6yuas/yr 2 / Vkpc/ 

Roughly speaking, this means that the astrometric survey can detect K = 1.1 times smaller 
accelerations than the RV survey and so can detect companions at K 1 / 2 greater separations 
or K 3 ^ greater periods. Naively, this would appear to mean that the fraction of grid stars 
with detectable accelerations (relative to the pre-RV sample) is given by 
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where Pf ~ P* and 

G(K) -> G (^) = log if (naive). (19) 

However, this simple treatment ignores the different ^-dependence in equations (|5|) and 
(14): the astrometric survey is sensitive to 2 components of acceleration, and corresponding 
to this, the mean value of sin 2 9 is twice as high as cos 2 6 over a sphere. In addition, the 
stars with low cos 9, which are most likely to evade the RV surveillance, are the most easily 
detected astrometrically. A more rigorous treatment yields, 

/accM = j dl °S r d ^ b gr min{y, yjl - (y/K) 2 }, y{r;m) = ( 20 ) 

Because the second factor differs significantly from zero over only a relatively narrow range, 
the first factor can be pulled out of the integral, which then yields equation ( |T8"D with, 

G{K) = J d\ogymm{y, ft - {y/Kf} = \og{K + + (21) 

and with P^ evaluated midway between the astrometric and RV sensitivities, 

d = k 3 / 8 P^ = . (22) 

'G(Mq + m) V2«thr/ 



Again, logPj- varies by only 0.6 over the whole mass range, so it is appropriate to evaluate 
/ace at (m) = 0.3 M Q . In this case, the correction for eccentric orbits is completely negligible 
because the effect is simply to slightly displace Pj. However, for the actual parameters of 
interest Pj- is near the peak of the DM91 distribution (where dfb/dlogP = 18/164 ~ 11%), 
so the evaluation of equation (|18D does not depend on the exact choice of Pj. For the fiducial 
parameters I have been using, / acc = 3.5%. Hence, a fraction / acc /(l — /rej,rv) ~ 5% of the 
grid stars that survive RV surveillance will have detectable astrometric accelerations. Here, 
I have used / re j irv = 30% to account both for eccentricity and for candidates eliminated by 
statistical fluctuations. 



4.2. Reference Stars 

Much of the foregoing derivation carries through for reference stars. There are three 
major differences. First, as I discuss in § 5.2, reference stars are likely to be closer: I adopt 
D = 600 pc. Second, the reference stars will be subject to measurements of much higher 
precision. Typically, target stars are expected to be observed relative to four reference 
stars at iV ~ 25 epochs, each time to a total positional precision cT pos ~ 1 //as. Since each 
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reference star is observed only 1/4 the time, and since their accelerations must be detected 
against each other, this implies that cr in equations (|T5|) and (116]) should be evaluated, 



VI /ias/ V25, 

a factor 6 times smaller than for grid stars. Consequently, K is a factor 6/0.6 times larger, 
K = 11. 

Third, depending on the ultimate design of SIM, the reference stars may be measured 
relative to the target only along one dimension ('ParaSIM') or in two dimensions ('Shared 
Baseline'). In the latter case, the analysis is identical to that given in § 4.1. Using 
K = 11 and equations (|T8| ) and (21), I then find / acc = 11%, implying that a fraction 



/acc/(l ~ /rej,rv) ~ 16% of reference stars that survive RV surveillance will have measurable 
accelerations. 

For the former case (which at the present time appears less likely), sin# in equation 



(14) should be replaced by sin#| cos0| where is a random angle on a circle. In principle 
this means that equation ( PD| ) should be replaced by a more complicated integration. In 
fact, this is unnecessary for the case K ^> 1, which is of relevance here. Since the RV and 
astrometric measurents are both one- dimensional, the fraction of stars with detectable RV 
accelerations at radius r rv will be exactly the same as those with detectable astrometric 
accelerations at r ast = i^ 1//2 r rv . Hence, if these two detection processes could be regarded as 
completely independent of one another, equations (|1~8|) and (p~9|) would hold. As discussed 



following equation ([19]), these detection processes are not generally independent. However, 
for K 3> 1, they are approximately independent because at the radii where the radial 
acceleration is detectable at all, there is only a very small chance that the astrometric 
acceleration will be undetectable. Applying equations (|TJ) and (JT9]), and taking K = 11, I 
find / acc = 9%, implying that / acc /(l - / r0 j,rv) ~ 12%. 



5. Implications 
5.1. Grid Stars 

In § 4.1, I showed that about 5% of SIM grid stars that survive RV selection will have 
measurable astrometric accelerations. If it were necessary to eliminate these stars from the 
grid well into the mission, then it would be necessary to build redundancy into the grid to 
"paper over" the resulting holes. In fact, this is not necessary. The fundamental reason 
is that, as I show below, while some stars surviving RV selection may have detectable 
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accelerations, almost none have detectable jerks. If this is the case, the accelerating stars 
can be fit to seven parameters (including two components of acceleration) instead of the 
usual five. From equation flT5|) this acceleration measurement decouples completely from 
the proper-motion measurement, and as discussed directly above equation (f[3), it decouples 
from the parallax measurement as well. Thus, the grid parallax and proper-motions (the 

main reasons for having a grid) are not significantly affected by fitting for acceleration. The 

1 li 

error in the mean position is increased by a factor cH /<To = 1.5. However, this error is more 
than an order of magnitude below any known requirement. For example, the positions are 
used for the tie-in to the radio reference frame, but the latter is known only to 20 //as, and 
the individual positions of quasars are much more poorly determined. Moreover, this slight 
degradation in positional error affects only 5% of the grid stars. 

What fraction of grid stars surviving RV selection will have detectable jerks? Noting 
that the parallax measurement again decouples (see § 4.1), I write the remaining positional 
dependence in one direction as ip(t) = ipo + Hot + (l/2)a ^ 2 + (1/ 6)jt 3 . I then evaluate bo- 
using equation (§) and invert the resulting (2 x 2) (proper-motionjerk) submatrix, thus 
finding 

*, = vwm% = ^(^)(^r. (24) 

yr 3 \2/ias/V5yr/ 
For circular orbits of period P, the astrometric jerk is given by 

._ m r /2vr\ 3 _ 8 /ias / m \/ r \~ 7 / 2 / £> V 1 
3 ~M^D\~p) yr 3 VO.lM JuOAU/ VkpcV ' ^ ' 

where to be conservative I have assumed that the jerk is in the plane of the sky and where 
I have written the formula in a way that is strictly appropriate only for m <C M . I 
will show below that this is the regime of the greatest concern. Hence, for the fiducial 
parameters, jerk is detectable at the 2 a level provided that r < 9AU(m/O.lM ) 2 / 7 . In 
the uniform-acceleration approximation, such accelerations escape RV detection only a 
fraction 14%(m/0.1 M Q )~ 3/7 of the time [see eqs. ©-(0)]. However, at r = 9AU, the 
uniform-acceleration approximation is already starting to break down and, as I discussed 
in § 3, for r < 4.5 AU, the RV survey almost never fails to detect stellar companions. I 
therefore estimate that for 1/2 dex in logP, the RV survey misses ~ 10% of grid candidates 
that go on to show marginally detectable jerk, i.e., about 0.4% of grid stars. 

To understand the effect of eccentricity, note that the worst case is a companion at 
periastron r of a highly eccentric orbit. Then the acceleration is the same as for a circular 
orbit at r, but the jerk is larger by \[2. Since stars spend little time at periastron and since 
the effect itself is small, I ignore it. 



Finally, I have ignored planetary companions of grid candidates. The plan is for the 
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grid to be composed of metal-poor K giants, which are not expected to have planets for 
several reasons. First, planet frequency among G stars (the progenitors of K giants) in 
the solar neighborhood is highly correlated with metallicity (Gonzalez 1997; Gonzalez, 
Wallerstein & Saar 1999). Second, planets occur in 47 Tuc at a rate that is much lower 
than that of solar-metallicity stars, and which is consistent with zero (Gilliland et al. 2000). 
Third, if gas giant planets grow from rock and ice cores, as most current theories suggest, it 
is difficult to see how they would get started in a metal poor environment. This optimistic 
assessment could prove wrong, but if so, it will become evident early in the RV survey of 
candidates, which would give plenty of time to modify strategy. For now, it is reasonable to 
suppose that such planets will be rare or non-existent. 

In brief, with the very reasonable fiducial parameters adopted here, I expect that 
RV screening can identify grid stars whose chance of being corrupted by an undetected 
companion is <C 1%. 

5.2. Reference Stars 

As I discussed in the introduction, contamination is a much more severe problem for 
reference stars than it is for the grid. First, the astrometric measurements are substantially 
more precise, meaning that they are more sensitive to contaminants. Second, the reference 
stars need to be brighter (to achieve this greater astrometric precision in a short exposure), 
which generally means that they need to be closer. This in turn increases the astrometric 
contamination, which scales oc D^ 1 . Third, their required surface density on the sky is 
higher, which means one cannot typically find metal-poor halo stars to serve as reference 
stars. Hence, a significant fraction is likely to have planets. Finally, uniform acceleration 
seriously undermines the function of reference stars whereas, as we saw in § 5.1, it is 
perfectly acceptable for the grid. Here I employ the formalism developed above to illuminate 
these problems. 

In order to be sensitive to very low levels of acceleration due to distant planets, the 
observer must find at least two reference starsQ that are themselves unaccelerating. This is 
because, if the target star and the reference star are found to be accelerating relative to one 
another, there is no way to determine which is truly accelerating and which (if either) is in 
uniform motion. It is only by finding two reference stars that are not accelerating relative 



1 In the case of ParaSIM, which measures relative offsets in only one direction, two reference stars must be 
found in each of two orientations. I ignore explicit consideration of ParaSIM here, but it is straight forward 
to extend the results presented to this case. 
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to one another that one can have reasonable confidence that it is they, and not the target, 
that are in uniform motion. If two such stars cannot be found, then one still has sensitivity 
to planets with periods P <^ T m , but not to outlying planets, including outlying companions 
of planets with short periods (i.e., planetary systems). 

The first point then, is that the RV survey must begin with enough candidates to 
have good prospects (say 95% probability) of finding at least two candidates that have no 
detectable RV acceleration. I showed in § 3, that with a 5 year, 60 ms" 1 survey, one could 
expect 30% rejection due to stellar companions and statistical fluctuactions. I have not 
evaluated planet contamination in this paper, but for illustration, I will assume that an 
additional 7% of candidates are eliminated by the RV survey due to planets. This means 
that if one wants to wind up with 2 reference stars, one should begin with 8 candidates. 

The second point is that if these are the only reference stars used, then the probability 
that one of them will prove to be an astrometric accelerator is high. In § 4.2, I showed that 
16% would have detectable accelerations due to stellar companions. However, there will be 
additional losses due to planets. For example, a Jupiter-mass planet at 3 AU would generate 
a velocity semi-amplitude of only 17sinims _1 , well below the detection threshold of the 
RV survey, but an astrometric semi-amplitude of 5 //as, quite easily detectable. Hence, an 
additional 5% loss due to planets is quite plausible. If so, the probability of at least one of 
the two reference stars being an accelerator would be 38%. 

There are basically four alternatives. 1) Accept that for more than a third of the 
target stars there will be sensitivity to closed planetary orbits, but not to distant planets. 

2) Accept fainter reference stars to increase D and so decrease astrometric contamination, 
but thereby degrade the astrometric precision and so the sensitivity to low-mass planets. 

3) Increase the number of reference stars to increase robustness. 4) Increase surveillance 
efforts to reduce contaminants. 

I have not much to say about the first two options, save that I would be disappointed 
if they were adopted. The third can be implemented at fairly low cost. If 8 candidates 
are initially surveyed, then there is an 87% chance that at least three will survive the RV 
survey, and 66% chance that four will survive. With three reference stars instead of two, 
the probability that at least two will survive rises from 62% to 89%. Of course, there is then 
also a high additional probability (39%) that one of these three will fai l, in which case the 
errors for detecting distant (but not close) planets would increase by y / 3/2. This is still a 
lot better than a total loss. The situation would be even more favorable with four reference 
stars. 



Finally, a more intensive RV study would be expensive in terms of big-telescope 
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time, but would be effective against both stellar and planetary companions. For example, 
reducing a TV from 60 to 20 m/s would decrease the number of stellar contaminants by 
A \riG(K) ~ 36%. It would not be possible to go below 20 m/s because this is the typical 
scale of photospheric fluctuations of K giants (Frink et al. 2001). However, if need be one 
could achieve the same effect with multiple measurements (assuming, as is almost certainly 
the case, that these fluctuations have power overwhelmingly on short timescales so that 
they do not couple to the acceleration measurements). 

6. Comparison with Monte Carlo 

Comparison with the Monte Carlo simulation of Frink et al. (2001) permits a direct 
check of the foregoing calculation, but first it is necessary to translate their parameterization 
into the one used here. They consider a RV survey with N=2 epochs separated by T rv = 5 yr, 
each with a vv = a/2 x 20 ms" 1 (including measurement errors and intrinsic instability of 
the stellar atmospheres). Hence, from equation (H), cr a = 8ms" 1 . I focus on the case 
of Xumit = ^ (their Fig. 6) corresponding to a tri r = VS&a- They demand that the rms 
scatter of the astrometric measurements be less than 1 /xas for a T m = 5 yr mission, which 
corresponds to a t hr = v / 720/ias/T^ l = 1.07/^asyr -2 , and place their grid stars at D = 2kpc 
(S. Frink, private communication 2001). From equation fll7l ) these values imply K = 2.23, 
hence G(K) = 0.67. From equation $2%), I find, logPj = 5.2, which I then apply to their 
adopted binary distribution function df b /d\ogP = 0.087exp[— (logP — 4.8) 2 /2 x 2.3 2 ], to 
find dfb/d\ogPj = 0.086. Equation ( |l8l) then predicts that 4.3% of their original sample 
should be found to be accelerators, or since 25% of that sample was rejected by their RV 
selection, 4.3%/0.75 = 5.8% of the grid stars. This compares with the value 5.5 ± 0.5% 
shown in their Figure 6, i.e., excellent agreement. 

7. Summary of Formulae 

I find that the fraction of grid-star or reference- star candidates that are eliminated by 
RV surveillance is equal to the fraction with binary companions having P < P*, where P* is 
given by equations ( |i~0"|) and ( |i~T| ) in terms of mean mass of the companion distribution and 
the threshold of acceleration detection, athr = Mcr a . Here <j a is the acceleration-measurement 
error, given by equation (|6|). (I used M = 2.) There is a small correction for eccentricity 
given explicitly by equation (|13"D, and of course some candidates will be falsely rejected due 
to statistical fluctuations, depending on the choice of M. 
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The fraction of initial candidates that survive RV surveillance but 
nevertheless have detectable astrometric accelerations is given by equation QIB"]), 
face — {3/4)(dfi,/dlogPj)G(K). The first term is simply the differential binary distribution 
evaluated at Pj, which latter is given explicitly by equation (|2"2"|). The second term is given 
by equation fl2~l"l) , G(K) = log (if + \/K 2 + 1) as a function of K, which characterizes the 
ratio of RV to astrometric sensitivies and is given by (|17|). In this case the correction for 
eccentricity is negligible. 

I find that a relatively modest RV survey can remove all but ~ 5% of grid-star 
candidates that will go on to show detectable accelerations, and all but ~ 0.4% of those 
that will show detectable jerks. I argue that it is only the latter very small fraction that 
must be eliminated from the grid. 

For reasons summarized in § 5.2, reference-star selection is much more demanding than 
grid-star selection. I find that to have a 95% chance of ultimately locating 2 reference stars 
(the minimum required to be sensitive to distant planetary companions of the target star) 
8 candidates must be initially surveyed for each target star. 
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pour L'Accueil et Les Echanges Internationaux. 
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